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Abstract. We consider some multivariate rational functions which have (or are conjectured to 
have) only positive coefficients in their series expansion. We consider an operator that preserves 
positivity of series coefficients, and apply the inverse of this operator to the rational functions. 
We obtain new rational functions which seem to have only positive coefficients, whose positivity 
, would imply positivity of the original series, and which, in a certain sense, cannot be improved 

C " ;) ■ any further. 
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1. Introduction 

Are all the coefRcients in the multivariate series expansion about the origin of 



1 — x — y — z — w+l (xy + xz + xw + yz + yw + zw) 

' positive? Nobody knows. For a similar rational function in three variables, Szego ^ has shown 

. positivity of the series coefficients using involved arguments. His dissatisfaction with the discrep- 

ancy between the simplicity of the statement and the sophistication of the methods he used in 
his proof has motivated further research about positivity of the series coefficients of multivariate 
rational functions. For several rational functions, including Szego's, there are now simple proofs 
for the positivity of their coefficients available. For others, including the one quoted above [ij, the 
positivity of their coefficients are long-standing and still open conjectures. 

In this paper, we consider the positivity problem in connection with the operator Tp {p > 0) 
, ^ ! defined as follows: 

OO . Tp: Rlxi, . . . ,a;„] ^ R|a;i, . . . ,a;„], 

(N ■ / 

r \ Vl^(l~P)2:i ' ■ ■ ■ ' l~(l-p)a;„ 



' By construction, the operator Tp preserves positive coefficients for any < p < 1, i.e., if a power 

, series / has positive coefficients, then the power series Tpf has positive coefficients as well, for any 

< p < 1. For example, via 

■ / 1 \ 1 

k>( , Tirtf - - , 

, — X — y — z + Axyz J l~x — y — z+j {xy + xz + yz)' 

positivity of the former rational function [2 implies positivity of Szego's rational function [7]. 
This is a fortunate relation, because the positivity of the former can be shown directly by a simple 
argument [1] while this is not as easily possible for the latter [5]. (Straub [5j gives a different 
positivity preserving operator also connecting these two functions.) 

This suggests applying the operator Tp "backwards" to a rational function / for which positivity 
of the coefficients is conjectured, in the hope that this leads to a rational function which again 
has positive coefficients, and for which positivity of the coefficients is easier to prove. We present 
some empirical results in this direction. Our results may or may not lead closer to rigorous proofs 
of some open problems. In either case, we also find them interesting in their own right. 
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2. Sharp Improvements 

Given a rational function /, we are interested in parameters p € [0, 1] such that T~^f has positive 
series coefficients. Because of = Ti/p, this is equivalent to asking for parameters p > 1 such 
that Tpf has positive series coefficients. Clearly, the set of all p > such that Tpf has positive 
coefficients forms an interval [0, Pmax) with a characteristic upper bound Pmax for each particular /. 
Computer experiments have led to the following empirical results. 

Empirical Result 1. Let f{x,y,z) ~ — x — y — z + Axyz). Let po be the real root of 
2x^ — 3x^ — 1 with pq w 1.68. Then po — Pmax- 



Evidence. (1) Pmax cannot be larger than po, because the particular coefhcient {xyz)Tpf — 
I + 3p^ — 2p^ fails to be positive for p> po- 

(2) CAD computations confirm that all terms {x'"'y™'z'')Tpf with < n, m, fc < 50 are positive 
for any < p < pg. 

(3) For p = 2430275/1448618, aU terms {x"y"'z'')Tpf with < n,m,k < 100 are positive. 
This p is the 15th convergent to po and only about 10"^^ smaller than this value. 

(4) For each specific choice of m,k, the terms {x"y™z'')Tpf are polynomials in n (and p) of 
degree m + k with respect to n. For < m, k < 10, CAD computations confirm that these 
are postive for all n > 1 and all < p < po- 

□ 

Empirical Result 2. Let f{x,y,z,w) = 1/(1 — x — y~z — w+ + xz + xw + yz + yw + zw)). 
Let Po be the real root of x'^ — 6x^ — 3 with po ~ 2.54. Then po — Pmax- 



Evidence. (1) Pmax cannot be larger than po, because the particular coefficient {xyzw)Tpf — 
3 + 6p^ — p"' fails to be positive for p> pa- 

(2) CAD computations confirm that all terms {x'^y"^ z^w'')Tpf with < n,m,k,l < 25 are 
positive for any < p < po. 

(3) For p = 730647/287378, aU terms {x"y"'z''w'^)Tpf with < n,m,kj < 240 are positive. 
This p is the 15th convergent to po and only about 10~^^ smaller than this value. 

(4) For each specific choice of m, fc, I, the terms {x"y™z'^w'')Tpf are polynomials in n (and p) 
of degree m + k + I with respect to n. For < to. A:, ^ < 5, CAD computations confirm 
that these polynomial are postive for all n > 1 and all < p < po. 

□ 



For the rational function / considered in Statement [21 our hope was dispelled that a direct proof for 
the positivity of Tp^^^^f could be found more easily than for / itself. However, some "suboptimal" 
values of p do lead to rational functions which have a promising shape. For instance, we found 
that 



1 



y — z — w + g (xy + xz + xw + yz + yw + zw) 
1 

1 — X — y — z — w + 2{xyz + xyw + xzw + yzw) + Axyzw 

Also note that it seems to be a coincidence that Pmax is determined by the coefficient of xyzw in 
Tpf, because this does not hold in the expansion of 

fix, y, z, w) = ■ ■ ■ 

l~x — y — z — w + ^[xyz + xyw + xzw + yzw) 

which is conjectured to have positive coefficients . Here we have Pmax < 1.66, by inspection of 
the coefficients {x"y"'z''w'^)Tpf for <n,m,k,l < 100, while {xyzw)Tpf = -§p* - §p^ + 6p'^ + 1 
is positive for p < 2.36. 
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Figure 1. 



3. Asymptotically Positive Coefficients 



Inspection of initial coefficients of 



1 — X — y — z — w + ^{xyz + xyw + xzw + yzw) 
max that become smaller and smaller as the amount of initial values taken 



■'real" value Pmax determined by the asymptotic behaviour of 



suggests values for p 
into consideration increases. Is the 
the coefficients for general pi 

Clearly, it is hard to extract conjectures about asyptotic behaviour by just looking at initial values. 
Instead, such information is better extracted from suitable recurrence equations by looking at 
the characteristic polynomial and the indicial equation of the recurrence jl]. Using computer 
algebra, obtaining recurrence equations for the coefficient sequences is an easy task. Often, the 
asymptotics can be rigorously determined from a recurrence up to a constant multiple fT, which 
cannot be determined exactly, but for which numeric approximations can be found. For instance, 
if p > Po (15 



\- 3v33)/2, we have 
(x"2/"^"«;")Tp(l/(l- 
(155 + 27%/33) 



X — y — z 



%j{xyz + xyw 



yzw))j 



K 



-4p + 3V33- 15)' 



3456 



in 



for K > 0.291 (Figure [T^ shows a„/ ((• • • )"n ^) for p = po + j^, supporting the estimate for K). 
This is oscillating. For 1 < p < pq, the asymptotic behaviour turns into 



k{i + Ip] 



3n 



(n —> oo), 

for K > 0.227 (Figure [TJd shows a„/ ((• • • )"n^^) for p = Po — jq, supporting the estimate for K). 
This is not oscillating, but ultimately positive. This supports the conjecture Pmax < Po ~ 16.1168, 
which is little news, however, as we already know Pmax < 1-66 by inspection of initial values. 
Other paths to infinity that we tried do not give sharper bounds on Pmax- So it seems that Pmax 
in this example is determined neither by the initial coefficients, nor by the coefficients at infinity, 
but somehow by the coefficients "in the middle" . 

We can consider asymptotic positivity of coefficients as an independent question which may also 
be asked for the rational functions considered in Statements [T] and [21 What are the values p > 
Pmax > 1 such that the series coefficients of Tpf are ultimately positive? Denote by p^ax the 
supremum of these parameters. We have carried out computer experiments in search for p^ax; 
and we obtained the following empirical results. 



Empirical Result 3. Let f{x,y,z) — 1/(1 



■ y — z + Axyz) . Then p° 



Evidence. Let e > (sufficiently small) and an^m.k '■= (a;"i/™z''')T2_e/. 

(1) First of all, we have p^ax — 2, because for e = 0, the asymptotics on the main diagonal is 

an,n,n if (-27)"n^2/3 ^) 

with K > 0.25, i.e., a„_„.„ is ultimately oscillating for e = 0. 
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(2) Let m, fc > be fixed and consider an,m,k as a sequence in n. A direct calculation shows 
that 



r=0 t=0 8=0 



n\ / n + m — r\ / n + m — 2r\ /n + m — r + t 
r I \ m — r ) \ s I \ t — s 



(2 - fA2(m+fc) 

= ^ + o(n'"+'=) (n ^ oo), 

m!fc! 

which is positive for n — > oo. 

(3) For arbitrary (symbolic) i > and the particular values < j < 3, the sequence a„^n+i,j 
satisfies a recurrence equation of order 3 which gives rise to 

a„,„+,,, ~i^(3-e)2"n-i/2 ^ oo) 

for some constants K depending on and e. Numeric computations suggest that these 
constants arc positive, and hence, an^n+i,j is positive for n oo. 

(4) For the particular values < «,j < 2, the sequence an,n+i,n+j satisfies a recurrence 
equation of order 3 which gives rise to 

a«,n+i,«+i ~ i^(3 - e)^"n"^ (n oo) 

for some constants K depending on i, j, and e. Numeric computations suggest that these 
constants are positive, and hence, a„^n+i.n+j is positive for n — > oo. 

□ 

In parts [3] and UJ we could not carry out the arguments for both i and j being generic. We did 
find a recurrence equation of order 6 for an^n+i.n+j for generic i, j, with polynomial coefficients of 
total degree 16 with respect to n, i, j, but this recurrence was way to big for further processing. 

5U1L '±. ucL j y jj ^ y ^ ,c ^ w ) — ^ I \ ^ — '■'^ — V — — w -r ' 

Then p^^^^ = 3. 



Empirical Result 4. Let f{x,y,z,w) — 1/(1 — x — y — z — w + + xz + xw + yz + yw + zw)) 



Evidence. Let p > I and an^m,k,i '■= {x'^y"^ z*'w^)Tpf . 

(1) First of all, we have pj^ax — 3 because for p > 3, the asymptotics on the main diagonal is 
determined by the two complex conjugated roots 



9 + 30p2 - 7pi ± 4p{p^ + 3)^6 -2p2 
9 ' 

Their modulus is {p^ — 1)^- As (p^ — 1)^ itself is not a characteristic root, it follows [S] 
that an,n,n,n IS ultimately oscillating for p > 3. 

(2) For fc > fixed, an,i,j,k is a polynomial in n of degree i + j + k + 1 whose leading 
coefficient is p2(i+j+fe) JHI j\lk\. Therefore an^ij^k is positive for rt ^ oo regardless 
of p. 

(3) For the particular values i = 0, 1 and < j,k < 2, the sequence an,n+i,j,k satisfies a 
recurrence equation of order 3 which gives rise to 

(p+V3)^ , . 

an,n+i,j,k ^ K — n-' 2 (n oo) 

for some constants K depending on i,j,k, and p. Numeric computations suggest that 
these constants are positive, and hence, an,n+i,j,k is positive for n ^ oo regardless of p. 

(4) For the particular values < j, fc < 1, the sequence an^n.n+j.k satisfies a recurrence 
equation of order 4 which gives rise to 

a„,„^n+j.k ^ Kil + p)'''"n"^ (n ^ oo) 

for some constants K depending on j, fc, and p. Numeric computations suggest that these 
constants are positive, and hence, an,n,n+j,k is positive for n ^ oo regardless of p. 
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Figure 2. 

(5) The main diagonal an.n,n.n satisfies a recurrence of order 4 which gives rise to 

a«,n,«,n (64 + _ 9) (2/ + 9p2 + 189 - 2(^2 + 3)^/2^)) "n-3/2 

for some constant K depending on p. 

Numeric computations suggest that K is positive for p < 3. For example, Figure [5] 
shows the quotients a„^„_„_„/((- • • )"ri^'^/2) fQj- p = S — jL., 

□ 



Stronger evidence in support of the conjectures made in the paper is currently beyond our com- 
putational and methodical capabilities. So are rigorous proofs. 
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